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Abstract 

We amplify previous arguments why mean curvature should be used as 
measure of integration in calculating the effective bending rigidity of fluid 
membranes subjected to a weak background curvature. The stiffening of the 
membrane by its fluctuations, recently derived for spherical shapes, is recov- 
ered for cylindrical curvature. Employing curvilinear coordinates, we then dis- 
cuss stiffening for arbitrary shapes, confirm that the elastic modulus of Gaus- 
sian curvature is not renormalized in the presence of fluctuations, and show 
for the first time that any spontaneous curvature also remains unchanged. 
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Introduction 



The fluid membranes of red blood cells and giant vesicles in water undergo mi- 
croscopically visible thermal shape fluctuations 0|. The thermal undulations lead 
to orientational decorrelation on a scale characterized by de Gennes' persistence 
length 1^. About fifteen years ago, one of us proposed that the decrease of cor- 
relation should be accompanied by a softening of the membrane. Specifically, the 
effective bending rigidity k' was predicted to obey 



kT 

k' = K- a — logM , (1) 
Stt 



where k is the bare bending rigidity and M the number of molecules making up the 
membrane. The numerical prefactor was found to be a = 1 in the first derivation 
and two later ones intended to support this result ||. In all three papers 
mean curvature was assumed to be the statistical measure. However, serious mis- 
takes were made in the calculations. Peliti and Leibler, whose result was the first 
to be pubfished |§, and all other authors [|[|, ||, |oHll|, |l| rederived eq. (0) for 
the renormalized bending rigidity, but with the prefactor a — 3. They took normal 
displacement of the membrane from its equilibrium shape as measure of integration. 
Only recently the problem was reconsidered [l^ when the mistakes were noted in 
the calculations employing mean curvature. In addition to correcting the previous 
calculation for the sphere [Q , a new calculation was done for the almost flat fluctu- 
ating membrane. It is based on a local bending mode expansion called hat model 
[Q. In both cases it was found that fluid membranes should be stiffened, not soft- 
ened, by their thermal undulations. The form of the relationship (|l|) is maintained, 
while the numerical prefactor changes sign, becoming a ~ —1. However, it has be 
mentioned that a few years ago Gompper and KroU found a membrane softening in 
Monte Carlo simulations of fluctuating vesicles p5| . 
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In the present work we deal once more with the problem of selecting the correct 
statistical measure, taking advantage of the analogies between polymers in two- 
dimensional space and fluid membranes. Subsequently, we consider a particularly 
transparent situation, cylindrical curvature, to rederive the effective bending rigid- 
ity in this case. Employing curvilinear coordinates, we then find the stiffening, i.e. 
the validity of eq. (|^) with a = — 1, for more general deformations of the initially 
flat membrane. In the same tensor notation, we confirm that the elastic modulus 
of Gaussian curvature is not affected by thermal undulations and show for the first 
time that this also holds for any spontaneous curvature of the membrane. 



Formulation of the problem 

The bending energy of a fluid membrane may be expressed by the Hamiltonian 

H = J (fcr (^^kJ^ - KC,J + kK^ , (2) 

where k and k are the bending moduli and Cg is the spontaneous curvature. The 
terms in parentheses represent the surface energy density and ^/g d^a is the local 
surface element of the membrane. Let X^a) , a = (cr^, cr^), be a parametrization of 
the membrane surface with the contravariant coordinates and cr^. Then d^a = 
da^da^ and g — det gij , where 

g,j = d,X ■ djX , with d,X = dXjda' , i , j = 1, 2 , 

are the coefficients of the first fundamental form. J and K denote trace and deter- 
minant of the extrinsic curvature tensor Jj — g^^Jkj: respectively, where 

J,j^N-d,jX, *,j = l,2 
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are the coefficients of the second fundamental form, N being (in the case of a closed 



surface) the outward unit normal vector to the surface. In other words, 
J = ci + C2 and K = ci ■ C2 , 

are twice the mean curvature and the Gaussian curvature, respectively, with ci and 
C2 being the principal curvatures. Note that with this convention J is negative for 
spheres and cylinders. 

The concept of an effective rigidity can be understood as follows: Starting from 
a basically flat, weakly fluctuating piece of membrane, we are interested in the free 
energy F needed to deform the piece into a fluctuating shape of nonvanishing back- 
ground curvature J. It turns out that the effect of the short wavelength fluctuations 
upon the elastic behaviour on large scales can be absorbed into the material pa- 
rameters of the bending Hamiltonian (|^). Accordingly, the contribution to F from 
the first term in (0) may be written as 



Here the curvature and the surface element refer to the nonfluctuating background 
and k' ^ K is the effective bending rigidity. 

As is common in shape calculations, the membrane is assumed to be unstretchable. 
The free energy F refers to a background surface of fixed area. The associated real 
membrane area is slightly larger because of the bending fluctuations. Any shrinkage 
of base area resulting from the fluctuations can be formally avoided by introducing 
a (fictitious) reservoir of basically flat membrane area. Alternatively, F can be taken 
to refer to the shrunken base area. Such a shrinkage does not affect the bending 
energy of the base if it is done by an isotropic scale transformation, as happens 
automatically in the case of a freely fluctuating spherical vesicle. Fortunately, these 
subtleties seem to be irrelevant in practice. 




(3) 
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The free energy of bending, T, can be obtained in two different, but equivalent ways. 
One may "freeze" the undulations in an instantaneous configuration and then cal- 
culate the bending energy of the rippled membrane at fixed undulation amplitudes. 
The squared fluctuation mode amplitudes are averaged only afterwards. 
Alternatively, one may start from the free energy of the fluctuating flat membrane 
and investigate how it changes when the membrane is subjected to a background 
curvature J ^ 0. The fluctuation dependent part of the bending energy F is now 
entropic, resulting from the effect of J on the mean square amplitudes of the undu- 
lation modes. 

In both cases we are faced with the problem of choosing the right measure of in- 
tegration. The undulations as expressed by the measure have to be kept constant 
when the membrane with "frozen" undulations is deformed in the large. Corre- 
spondingly, the entropic part of the free energy of background deformations has to 
be calculated from the change of the mode mean square amplitudes in terms of the 
measure. We will expand now on some of the arguments why mean curvature, or 
any quantity proportional to it, is a correct measure of integration. The scale of 
the measure is of no interest because we are dealing with differences, not absolute 
values of entropy. 



Considerations on the measure of integration 

A heuristic argument for choosing mean curvature as statistical measure is the fact 
that it can be regarded as the physical strain of the fluid membrane, much like 
curvature in the case of a (stiff) polymer. Accordingly, the superposition of two 
deformations, e.g. uniform background curvature and a (sinusoidal) ripple, begins 
with adding the strains as they are in the absence of the other deformation. The 
simultaneous presence of both deformations will produce coupling terms in the ef- 
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fective total bending energy that are quadratic in both the background curvature 
and the ripple amplitude. Employing the strain as statistical measure has the par- 
ticular attraction of avoiding spurious coupling terms. They arise from a change of 
the ripple amplitude in terms of the strain when other measures are used such as 
normal displacement from the background or reorientation of the layer normal. 

From a technical point of view, mean curvature can serve as measure of integration 
only if its local values determine the shape of the polymer or membrane, respec- 
tively. A one-to-one relationship between shape and local curvatures (which may 
be molecular) clearly exists for the polymer in the plane, a simplified model system 
frequently considered in the following. For membranes, a similar invertiblc map- 
ping is particularly transparent in the hat model which expands the fluctuations of 
the infinite planar membrane into local mean curvature modes p3[. There are only 
minor restrictions to the local mean curvature fiuctuations in other geometries, 
such as periodic boundary conditions, cylinders or spheres (see below). However, 
any (small) membrane deformation compatible with the boundary condition can be 
completely specified either by local normal displacements or mean curvatures. In 
all these cases Gaussian curvature fiuctuates together with mean curvature, but the 
former is completely slaved by the latter because of the boundary conditions. The 
general possibility of invertiblc mappings between shape and mean curvature, in- 
cluding local area conservation and lateral flow, has been discussed elsewhere p^ . 

Finally, there is a cogent reason to select (mean) curvature as statistical measure 
for polymers in 2D or fluid membranes. It stems from the fact that in both systems 
every local strain fluctuates independently of all the others (apart from the minor 
restrictions just mentioned). As a result, the total fluctuation entropy is the sum of 
the equal entropies of all the molecular (or other local) fluctuations. Independence 
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and equality imply that the entropy Sq of a sinusoidal (mean) curvature mode on a 
flat background does not depend on its wave vector q. This is because for all g 7^ 
the molecules are in the same way homogeneously distributed over the phases of 
the sine wave. The g-independence of Sq indicates that mean curvature is a suit- 
able measure of integration since the mean square curvature amplitudes also do not 
depend on q. Let us add as a remark that the absence of a wave vector dependent 
Jacobian can perhaps be regarded as characteristic of a correct measure. 
The condition of the mode entropy being independent of q is violated if the entropy 
is calculated in terms of other statistical measures. Taking normal displacement, 
one generally derives 

Sq = -fclogq^ + const. , 

from the variation with 1/q* of the mean square displacement amplitude of a si- 
nusoidal fluctuation mode. The g-dependent part of Sq is half as strong if normal 
displacement is replaced by the tilt angle of of the polymer or membrane with re- 
spect to the background. Evidently, neither normal displacement nor tilt angle 
is a correct measure of integration. It should be noted that in many statistical- 
mechanical calculations the choice of the measure is irrelevant. However, it matters 
whenever the coupling of two bending deformations is considered as will be demon- 
strated below. 

After these general statements let us address particular model systems to deal with 
the statistical mechanics of both polymers in a plane and fluid membranes. They 
will conflrm the claim that mean curvature is a correct measure of integration for 
fluid membranes. Examples of the invertible mapping between membrane shape 
and local mean curvature will also be given, including the restrictions on the latter. 



We first have a look at the comparatively simple problem of a polymer chain fluctu- 
ating in a plane [0 . The configurational space of the polymer is swept completely 
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by varying the angles between adjacent links. Disregarding fluctuations of link 
length, one may speak of local curvatures instead of angles. Each joint between 
subsequent links has its own free energy arising from its curvature fluctuations. Al- 
though the change of a particular angle rotates the two sides of the polymer chain 
with respect to each other, it does not affect their free energies. We assume here 
that the polymer is almost straight. In a continuum description such a polymer is 
equivalent to a (thin) rod or wire. 

Using local curvature as the measure of integration in dealing with fluctuating poly- 
mer chains is not in conflict with the fundamental measure which for mass points 
is translation. Local curvature is a shorthand notation for relative displacement if, 
e.g., each joint of the polymer chain is regarded as a mass point. Then a (small) 
normal displacement of point n + 1 from the tangent to the polymer defined by 
the two preceding points is proportional to the curvature at point n. In contrast, 
the absolute normal displacement of a nearly straight fluctuating polymer from a 
base line is not a suitable statistical measure. It implies a freedom of motion for 
each individual mass point which actually is suppressed by cohesion and bending 
stiffness. The remaining freedom is fully accounted for by the curvature. 
The polymer allows a simple demonstration of the effect of a wrong statistical mea- 
sure on the curvature amplitude of a ripple. We consider a sinusoidal wave on a 
polymer with a uniform background curvature — 1/r, the negative sign indicating a 
downward bend. A sketch of the situation is given in Fig. 1. Let us keep all normal 
displacements i/(s) fixed, s being the arc length of the background, while we bend 
the polymer as a whole. This amounts to using the displacement measure. The 
procedure results in a weakening of the ripple curvature because the total curva- 
ture changes from d^u/ds'^ to —1/r -|- d^u/ds'^ + u/r^ up to linear order in u, with 
the two terms constituting the ripple curvature being of opposite sign. Clearly, the 
change of the ripple curvature amplitude is physically not acceptable so that normal 
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displacement does not appear to be the correct statistical measure. 

Continuing the two-dimensional polymer into the third space dimension results 
in a bent membrane with zero curvature in the new direction. The similarity of 
such a membrane with the polymer in a plane strongly suggests that some kind of 
curvature rather than displacement from a fixed base should be the right measure 
of integration for fluid membranes. 

In order to extend the preceding argument to membranes that can be bent not 
only in a single direction let us consider a fluid membrane made of mass points. 
A particularly simple situation arises for the almost flat, weakly fluctuating mem- 
brane: Like the polymer, the idealized fluid membrane can be build up by adding 
the constituent mass points one by one. For easy statistics, we may think of a (lo- 
cal) lattice of mass points, quadratic or hexagonal, which is formed row after row. 
A new molecule continuing a straight sequence completes a plaquette comprising a 
certain number of preexisting neighbours. If the number is large enough the pla- 
quette can be used to define the mean curvature at its center. Five molecules is the 
minimum; they may be thought to form a centered square, as is illustrated in Fig. 2. 
The height of the last molecule added, relative to a suitably chosen tangent plane, 
finalizes the value of the curvature J. On a quadratic lattice, J may be taken (in 
lowest order) to be the sum of the (normal) curvatures of the two diagonals of the 
plaquette. Let us emphasize, that it is not the relative height of the final molecule 
that gives a parametrization of the phase space of the particle in the middle of the 
plaquette, but the simultaneously varying mean curvature. This quantity is invari- 
ant with respect to the direction along which rows are formed one by one. The 
tangent plane may be replaced by the base plane, shifted to the local height, if the 
fluctuations are weak enough to leave the membrane almost flat. Obviously, the 
last simplification is possible in the example of Fig. 2. 
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Regardless of the measure of integration, it is convenient to describe the thermal 
undulations superimposed on a background shape of given J and K in terms of 
normal displacement i^((t). To linear order, v and the associated extra curvature 
— J' — J are related through 

i^[J^-2K + D^] V , (4) 

being the Laplace-Beltrami operator (see eq. below and subsequent defini- 
tions). The difference ^ is equivalent to J as statistical measure. 
Undulation modes can be defined as eigenfunctions of the operator represented by 
the bracketed expression. For surfaces of uniform background curvature, such as 
planes, cylinders and spheres, this ensures that for each mode the functions ^((i) 
and v(ci) are proportional to each other in a first order approximation, the factor 
of proportionality depending on the wave number. The transformation can be 
inverted whenever for a given undulation mode the operator in brackets does not 
vanish or affect the membrane area linearly in the mode amplitude. Well-known ex- 
amples of such exceptions are the three I — 1 modes of the fiuctuating sphere which 
actually represent translations and the I ~ mode which requires area changes 
of the unstretchable membrane. With periodic boundary conditions the only lost 
mode is uniform mean curvature (q = 0) as it is incompatible with the bound- 
ary conditions. In the case of the cylinder two long-wavelength curvature modes 
are lost because they represent translations normal to the cylinder axis. Periodic 
boundary conditions along the axis and area conservation suppress two other mean 
curvature modes. The loss of modes implies that the local mean curvatures are not 
totally free to fluctuate independently of each other. However, these restrictions 
seem irrelevant in the present context for membranes comprising large numbers of 
molecules. 
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Basic formulas 

In order to calculate the bending Hamiltonian (|^) for a fluctuating membrane, 
we first present expressions for the mean curvatures and the surface elements of 
membrane shapes lying close to a given background shape. The latter is also referred 
to as the projected shape. 

If X{(f) ,(7 — (ct^jCt^), is a parametrization of the background configuration, any 
slightly deviating configuration X'{d-) may be represented in normal gauge |Q, |l^, |l^ 
by 

X'ia)=X{a) + iy{a)N{a) , (5) 

7V(ct) being the local outward unit normal vector to X{d-). From ^ we find that 
up to quadratic order in v the surface element of the deformed membrane obeys 
(cf.@,|[000) 

^/^^V9{i-'^J+l{Diy)^+Ki.^) . (6) 

In the same approximation the curvature J' becomes 

f ^ J + [D^ + -2K) f J [Dvf 

+,f^ D.vDjV + (J3 - 3Jif ) + 2J'^ vDiDjV + 

+vDr,vg^^{D,J^^ + JTl), (7) 

Here, gij and are the coefficients of the first and second fundamental forms 
of the background surface, respectively. are the Gauss-Christoffel symbols of 
the second kind derived from gij. Di is the covariant derivative being defined as 
DiV = diV — dv/dcr^ , and DiDjV — dijV — T^jdki^. = g^^ DiDj is the Laplace- 
Bcltrami operator and we use the abbreviation {Dv)^ = DV Dii'. The variation of 
the Gaussian curvature K need not to be considered in (|^) thanks to the Gauss- 
Bonnet theorem 
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Renormalization of k for the cylinder via frozen 
mode amplitudes 

Let us now derive the effective bending rigidity for cylindrical background curvature 
in the "frozen mode" approach. We first write down reduced versions of (^) and 
that apply to this special case. Using geodesic coordinates x = r ■ (p and z to 
parametrize the cylinder, where ip and z are usual polar coordinates, we find with 
J = -l/r: 



/^=l + lu + l{Vuf +0(1^'') (8) 
r 2 



and 
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r r 2 r 

l.„ .1 „ 1 



^ {d^iyf - 2-iyd^.,iy - -^ly^ + 0(z.3) . (9) 

The linear part of J' in (||) leads to the transformation and its inversion for the 
cylinder, 

(l/r2 + A)j/ = C ^ = A-\l + l/r^A-^)-^C , (10) 

where A is the usual Laplacian for a flat base. Obviously, the thermal undulations 
on a cylinder can be expanded, like those on a flat background, into sinusoidal waves 
which at the same time are curvature and normal displacement modes. 
We will consider a quadratic piece of membrane subjected to a weak cylindrical 
background bend and satisfying periodic boundary conditions. Alternatively, the 
membrane may be closed to form a cylinder so that the edge of the square equals 
the circumference in length. In the displacement representation the sine and cosine 
deformation modes read 

v — a^siii{q f) and v — bq'COs{q f) , (11) 
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where r — {x. z) and q = (q^, qz) with qx = 27:71/ L and qz — 2Trm/L. The numbers 
n and m are running from ~L/{2a) to L/ {2a), but only over half the q plane. The 
constant mode (n, m) — (0, 0) is excluded. L denotes the size of the background 
shape and a is the molecular distance in the membrane. The operator acting on v 
in ([l0| ) is ~ 9 ^) for a given mode. 

In order to renormalize the bending rigidity let us first consider the effect of a 
single sinusoidal mode. A straightforward calculation starting from (^) and and 
involving a partial integration then results in the Hamiltonian 



1 ,2 ^ 

— kJ V q'da;dz 

2 ^ ^ 



1 



-'- I l/'i / 2 ^2\2 2,11,' -'^2,32 -'^\2 



base 



(12) 



where Abase = denotes the area of the nonfluctuating background shape. Terms 
linear in v do not occur as they vanish on integration. Of the three terms in the 
brackets, the second one equals the mean value of thus giving the regular bending 
energy of the mode. Note that if the curvature measure applies, the curvature 
difference ^ is kept constant in the " frozen mode" approach when the membrane is 
deformed. Accordingly, the third term is due to the coupling of the mode curvatures 
to cylindrical curvature. 

The statistical averages of the squared mode amplitudes are obtained from the 
equipartition theorem, 



Inserting (|T^) in (O) and absorbing the third term of (O) in the first leads to 



2 / 2 \ [l/r^-q-^y AbaseK J H 2 
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Adding up the contributions of all fluctuation modes to the renormalization of k 
yields to lowest order in kT/n 

where is omitted because oi q > 1/r. The sum is assumed to be small as 
compared to unity. It does not contain (n, m) — (0, 0) and some seemingly divergent 
low wave vector modes. The upper limits of n and m are large numbers. Replacing 
the sum by the double integral 



we arrive with good accuracy at 



kT 

k' = kH logM , (16) 

Stt 

M being the number of molecules in the membrane (or half this number for the 
bilayer). Evidently, we have recovered the result derived previously for the fluctu- 
ating spherical vesicle, i.e. eq. (|l]) with a — —\. 

The effect of thermal undulations on the bending rigidity of fluid membranes was 
attributed to a coupling between the mode curvatures and the background cur- 
vature. To confirm the role of coupling, we may for a moment compensate the 
uniform background curvature by a spontaneous curvature. For the cylinder this 
means Cs = —1/r and J' — Cg = ^. Accordingly, the surface energy density is sim- 
ply \'^^ — 5^ ((77 + ^) up to quadratic order in ^. The additional quadratic 
terms arising in the case of interest, Cs = 0, can therefore be regarded as coupling 
terms that give rise to a renormalization of the bending rigidity. 
For a direct understanding of the coupling let us consider two special types of rip- 
ples, those exactly parallel or orthogonal to the cylindrical background curvature. 
Parallel ripple: When the ripple is parallel we can use the analogy to a polymer 
in a plane. This makes it attractive to use a more direct approach than the general 
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calculation based on ^{x,!/). Starting from a straight and stiff polymer or a wire, 
we impose on it a periodic deformation 

£,{s) = Cg sm{q s) , 

where ^ is curvature, s arc length and q wave vector. The undulation diminishes 
the projected or background length of the wire. Alternatively, it increases the total 
arc length L' of the wire (when connected to a reservoir) at fixed projected length 
L. In both cases the ratio of the two lengths is 

L'/L = 1 + ii^ (17) 

to lowest order in i^. 

A frozen bending mode of the polymer is similar to a wire with a permanent wave. 
Any sinusoidal bending deformation allows a given angle difference between the ends 
of the wire to be spread over L' instead of L. This reduces the associated bending 
energy density by {L/L')'^ and the total bending energy by L/L' as compared to 
a polymer or wire without waviness. The analogous ratio for parallel ripples in 
membranes is A/A' , where A' and A are the real and projected areas, respectively. 
Absorbing the effect into a renormalized bending rigidity, we have k' = k A/ A'. 
This agrees exactly with what we obtain from (^4|) for a single mode of type {qx, 0) 
if we make use of A' = A{1 + jq'^ a^^ q). The negative sign in front of the correction 
term in k' indicates the softening in this case. 

An amusing illustration of the effect of parallel undulations is a rather stiff polymer 
closed to form a ring. Its thermal undulations in terms of curvature are not affected 
by the uniform bend and vice versa. Only the decrease of the ring radius in the 
presence of undulations signals a lowering of the effective bending stiffness of the 
polymer. 

Orthogonal ripple: Ripples orthogonal to the cylindrical base curvature stiffen 
the membrane, in contrast to parallel ones. This is due to the fact that at the 
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crests of such a ripple the extra curvature ^ is of the same sign as the cyhndrical 
curvature, and at the same time the membrane area is increased as compared to 
the background shape. In hnear approximations the two effects obey 

J.Z r J 

and 

r 

Calculating from these two expressions the total bending energy of a membrane 
with a single excited mode of wave vector (0, g^), we find the renormalized bending 
rigidity of the cylinder to be 

K' = nil-i-ql)alJ. (18) 

According to this formula the orthogonal ripple stiffens the membrane. It is four 
times as effective as a softening, parallel ripple of equal wavelength. 
There is a second correction quadratic in i/ to the bending energy due to the coupling 
between background curvature and orthogonal ripple. It originates from the sloped 
regions of the ripple where the membrane partially escapes cylindrical curvature. 
The reduced principal curvature is {cos(p)/r — (l/r)(l — + 0{lp^)), ip being 
the slope angle. For a single mode of type (0, qz) the second effect and the increase 
of membrane area due to the slope result in a softening given by 

«' = '^(l - «o,,J • (19) 

Combining the opposite rigidity corrections of (^8|) and (^9|) results in a stiffening 
which agrees with Q. 

At this point let us briefly consider what would be different in a derivation of 
k' based on the displacement measure. Keeping v fixed instead of ^ while the 
background is cylindrically bent implies that the regular part of the bending energy 
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density is 1/2 k(Az^)^ instead of 1/2 k{/S.v — v/r^Y'. Accordingly, the single mode 
Hamiltonian would have to be recast in the form 

/— nJ'^ \fq' Axdz = 
2 

if again the third term is to represent the coupling of the undulations to cylindrical 
background curvature. The negative signs before and in the coupling term 
suggest that both parallel and perpendicular ripples soften the membrane. Absorb- 
ing the third term into the first renormalizes k. The same manipulations as above, 
including thermal averaging and summation over all modes, would lead to 

kT 

k' = K-3 — logM . 

This is the result obtained by all other authors. They generally used the displace- 
ment measure and took the approach via undulation mode entropies. 
The undulations of unstretchable fluid membranes must be accompanied by lat- 
eral displacements. They can be analyzed mode by mode and consist of local and 
global parts. The local part redistributes membrane material between the crests 
and troughs of a ripple in a bent surface and in general varies linearly with the 
mode amplitude. The global part takes account of the area absorbed by the un- 
dulation and goes with the square of the amplitude. We have disregarded lateral 
motion, with one exception: In the case of the wavy polymer it is automatically 
included by employing arc length as independent variable. Inspection shows that 
lateral displacement need not be taken into account in calculations of the effective 
bending rigidity of polymers and membranes in the usual approximation that is 
quadratic in the mode amplitudes . 
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Renormalization of and k via undulation mode 
entropies 

The effective bending rigidity has been derived above for uniform cyUndrical cur- 
vature and earUer for the spherical membrane, both times in the "frozen mode" 
approach. The same stiffening was obtained in both calculations. Let us try now to 
extend the result to more general geometries, showing simultaneously that the un- 
dulations do not affect the spontaneous curvature Cs and the modulus of Gaussian 
curvature, of the membrane. Obviously, we have to limit ourselves to "com- 
pact" shapes such as squares, circles and spheres. For the following calculations it 
is advantageous to follow the approach based on fluctuation mode entropies. The 
curvature increment ^ = J' — J will again be used as measure of integration. 
In order to calculate the free energy of the fluctuating membrane and to renor- 
malize its elastic parameters, we employ a field-theoretic approach as has been 
introduced by Forster in this frame. Considering some background configura- 
tion with parametrization X{a), which describes the effective shape, we expand the 
bending Hamiltonian in (|^) to second order around X{a). We drop linear terms, 
assuming the background to be an equilibrium shape. The free energy F of X may 
then be written as 0, 

Fix) = H{X) + IfcTlogdet J^"^^^ + 0{{kTf) , (20) 

where ^ parametrizes the configurations X' in the vicinity of X. After expanding 
the r.h.s. in powers of background curvature one may introduce effective elastic pa- 
rameters marked by a prime, such as k', and an effective Hamiltonian H' employing 
the new parameters in conjunction with J and K of the background shape X. 
To perform this program we first have to calculate the bending energies of the 
fiuctuating membrane. Starting from eq. we obtain 

J'2 = j2 + 2 J(L>2 + - 2K)v + ((1)2 + j2 - 2K)vf 

18 



-J^ {Dvf + 2J.r^ D.vDjV + 2J(j3 - ■iJK) v"^ 

vD^DjV + 2Jv Dnv g'HD^Jl' + J T^) + 0{iy^) . (21) 

Combining this with (^) yields 

-2j'^v{D^ + J^- 2K)v 

+AJJ^'vD,DjV + 2jrW,iyDjiy + {{D^ + - 2K)vf 

+JKiy^ + 2J(J3 - '3JK)iy^ + 2JiyDniyg'^ {DjJl' + JT'^^) + 0{iy^)]y/^ . 

(22) 

We also need 

J'y^ = [J + {D^ + J2 - 2K)i^ ~ 

-Jv{D^ + - 2K)v + 2.r^vD^i:)jV + .r^D.vBjV 

+JKi^^ + (J3 - 3JKy + uDr,ug''{DjJl' + JF^^-) + Oiv"^)]^ . (23) 

In the following the background curvatures J and K are assumed to be weak and 
slowly varying. For all wave vectors except for the smallest ones the predominant 
terms quadratic in v are those that contain at least two derivatives acting on z/. In 
calculating the effect of the short wavelength fluctuations upon the long distance 
elastic behaviour of membranes comprising a large number of molecules, it is suffi- 
cient to take into account only the latter. Giving rise to cutoff dependent effects, 
they are the only ones that can contribute to the renormalization of the elastic 
parameters. To replace the normal displacement by ^ we need the inversion of (^ 

1^ ^ D-^{1 + {J^ -2K)D-^)-^^ , (24) 

In the spirit of the above simplification, ( p^ can be approximated for terms which 
are quadratic in ly at the stage of J' and g'^^^ (see eqs. (||) and (0)) by 
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Disregarding irrelevant terms, we may thus write the bending energy as 



-KCs • (J + C + ^{-JD-^ + rWiDj/D'^)0 + irrel. terms] 



(25) 



where we have performed some partial integrations, assuming either periodic bound- 
ary conditions or undulations limited to a patch of membrane. 
We are now in a position to set up the functional integral sweeping curvature space 
that gives the effective free energy F of the background shape X with fields J and 
K. In mechanical equilibrium linear terms in ( p5|) drop out so that we may write 



F 



fcTlog J PKlexp J d2a^(|j2-Kc,J 



D- 



2JJ'^D.Di 



2J2 



2c, 



2 D4 



D4 



irrel. terms I S, 



(26) 



Extracting the background energy and integrating the multiple Gaussian we are 
led to 



F 



d 

kT 



(f 

Tr' log 



J2 



2kT 



2D^ 



2J.rW,Dj 



2J2 



2c. 



(27) 



where / is the identical operator. The prime on Tr indicates exclusion of some 
modes at the low wave number cutoff. We assume rather weak background bends, 
such that 



(28) 
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and may thus expand the logarithm in powers of them. A further hmitation of 
may arise from Cs ■ The traces occuring in ( p7| ) are known to be [ p^ p3[ : 



[ d^a Jg ■ — log - and (29) 
J in a 

|dVV5^^~log^, (30) 



L and a being upper and lower cutoffs due to finite membrane size and molecular 
distance, respectively. They are easily calculated for uniformly bent shapes such as 
squares, spheres and cylinders. One obtains 

-|dVVgc..J.(l + |^log^) . 



(31) 



The middle term of ( ^l| ) is equal to the free energy of the corresponding flat reference 
membrane. After subtracting it we are left with the free energy of bending which 
may be written in the renormalized form 

r - 1 dV {]^n'J^ - k'c'J) , (32) 



where 



/ kT , L , , 

K ^ K+— log- 33 
47r a 

c's = Cs . (34) 



For comparison with ([l]), note that L'^ /a? = M. Evidently, the effective bending 
rigidity comes out exactly as calculated for the sphere ||l^ and the cylinder, while 
the spontaneous curvature is not modified by the thermal undulations. 
Formula (^2|) omits the Gaussian curvature term of the bending energy. We dropped 
it in the unrenormalized bending energy since, as a consequence of the Gauss- 
Bonnet theorem, the variation of Gaussian curvature does not affect its integral if 
the membrane is closed or satisfies periodic boundary conditions. Nevertheless, an 
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additional Gaussian curvature term could have emerged in the calculation of the 
effective bending rigidity, as it does when normal displacement is the measure of 
integration . From the fact that such a term is absent in eq. (^l|) we may infer 
that R is not affected by thermal undulations, i.e. 

R' = R , (35) 

in agreement with previous considerations |l3| . 

In the preceding calculation of the effective bending rigidity we did not specify the 
shape of the membrane, except for stipulating that it be "compact". Although the 
formula for k' seems to be rather robust, because of the logarithmic dependence on 
wave vector cutoffs, it is not expected to hold for filamentous membranes. Specifi- 
cally, the bending rigidity of the membrane closed to form a cylinder will not stiffen 
indefinitely as the cylinder is made longer. This is why we chose periodic boundary 
conditions with an axial period as long as the circumference in our calculation for 
the cylinder. 

The renormalized elastic parameters calculated above are independent of the base 
curvatures J and K. It is an open question to which extent they remain valid if 
these curvatures are nonuniform over the base area. In the case of a wavy back- 
ground one commonly uses as upper limit in its wavelength or a fraction thereof 
instead of the membrane size L. This is because ripples of larger wavelength are 
part of the immutable background curvature. 

Conclusion 

We have renormalized the bending elastic parameters of fluctuating fluid membranes 
for cylindrical and more general membrane shapes. The same logarithmic increase 
of the bending rigidity with membrane size was obtained as previously for spherical 
vesicles. In addition, the modulus of Gaussian curvature and, for the flrst time, 
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spontaneous curvature were shown not to be affected to lowest order by thermal 
undulations. Apart from the illustrative treatment of the cylindrical membrane, we 
used arbitrary curvilinear coordinates in the calculations to obtain generally valid 
results as far as possible. 

The predicted stiffening of fluid membranes by their thermal undulations may seem 
contrary to intuition, although there is no conflict with de Gennes' |^ derivation of 
the persistence length of membrane orientation. In the case of polymers the loss of 
orientational correlation unavoidably results in a breakdown of bending elasticity. 
This could be different in the case of membranes because of their two-dimensional 
connectedness and, in particular, the well-known scale invariance of bending en- 
ergies. Moreover, it remains to be checked if the stiffening survives at very low 
bending rigidities (k ~ kT), where the nearly flat approximation is inapplicable at 
any scale. 
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Figures 




Figure 1: Wavy polymer with fixed ripple curvature amplitude on flat and bent 
background, respectively. Note the enhanced normal displacement from the back- 
ground on the right-hand side. 




Figure 2: Lattice model of a membrane fluctuating around a flat base. The heavy 
lines are grid lines and shown only up to the membrane edge. The thin lines indicate 
membrane shape or, at membrane edge, normal displacement from the base. The 
five square-shaped points determine the mean curvature in the middle (solid square) 
of the plaquette 
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